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1 Introduction 

The subject of this paper is a special class of configurations, or patterns, of 
intersecting circles in constant curvature surfaces. The combinatorial aspect of 
such a pattern is described by a cellular decomposition of the surface. The faces 
of the cellular decomposition correspond to circles and the vertices correspond 
to points where circles intersect. (See figures g and §.) 

In the most general case that we consider, the surface may have cone-like 
singularities in the centers of the circles and in the points of intersection. 

In particular, we treat the problem of constructing such circle patterns when 
the cellular decomposition is given combinatorially, the intersection angles of the 
circles are prescribed by a given function on the edges (and, possibly, cone an- 
gles are prescribed by functions on the faces and vertices). Using variational 
methods, and the 'method of coherent angles,' we prove existence and unique- 
ness results. The most fundamental one of these is theorem 3L Similar results 



were o btained by Bowditch | Bow91 |, Garrett | Gar92 |, Rivin Riv99 |, and Leibon 
[ LeiOl ] . Our theorem treats circle patterns in surfaces which may be flat or hy- 
perbolic, closed surfaces or surfaces with boundary, and there may be cone-like 
singularities in the centers of circles and in their points of intersection. From it 
we deduce Rivin's theorem on ideal hyperbolic polyhedra [ Riv96|| , and theorem 
[|, which is the analogous theorem for higher genus surfaces. This higher genus 
analogue is new. 

Our method is based on two new functionals — one for the Euclidean and 
one for the hyper bolic ca se. We sho w how t he functionals of Colin de Verdiere 
[CdV91|, B raggcr | Bra92 |, and Rivin [Riy94] can be derived from ours. Leibon's 
functional |Lci01| seems to be related as well. 



2 A brief survey of relevant results 

A circle packing is a configuration of disjoint discs which may touch but not 
intersect. In 1 936, Ko ebe published the following theorem about circle packings 



in the sphere Koe36 | 
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Figure 1: A cellular decomposition of the torus and the corresponding circle 
pattern with orthogonally intersecting circles. 
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Theorem (Koebe). For every triangulation of the sphere there is a packing 
of circles in the sphere such that circles correspond to vertices and two circles 
touch if and only if the corresponding vertices are adjacent. This circle pattern 
is unique up to Mobius transformations of the sphere. 

Any graph without loops or double edges, embedded in the sphere, is con- 
tained as subgraph in the 1-skeleton of some triangulation. Hence, any such 
graph is the touching graph of a circle pattern in the sphere. Of course, the 
circle pattern is in general not unique up to Mobius transformations. 

There is, however, a way to generalize Koebe's theorem to cellular decom- 
positions of the sphere whose faces are not necessarily triangular without giving 
up the uniqueness part. The right class of cellular decompositions to consider is 
the class of cellular decompositions that are combinatorially equivalent to con- 



vex polyhedra. By a theorem of Steinitz Ste22 , SR34], these polytopal cellular 
decompositions are characterized by the following two conditions. First, the cel- 
lular decomposition is supposed to be regular, which means that that closure of 
each open cell is a closed disc. For the 1-cells (edges) this means that loops are 
not allowed. For the 2-cells this means that there are no identifications on the 
boundary. The second condition is that two closed cells intersect in one closed 
cell, if at all. This forbids double edges. More generally, it implies Steinitz' 
Bedingung des Nichiibergreifens: If each of two faces is incident with each of 
two vertices, then there is an edge which is incident with the two faces and the 
two vertices. 

So given a polytopal cellular decomposition of the sphere, there is a cor- 
responding circle packing. Uniqueness is retained if an additional condition is 
imposed on the packing: Consider the vertices incident with a face of the cellu- 
lar decomposition. The corresponding circles touch in points which correspond 
to the edges of the face. It is required that these points lie on a circle which 
intersects the original circles orthogonally. This is the content of the following 
theorem. 

Theorem 1. For every polytopal cellular decomposition of the sphere, there ex- 
ists a pattern of circles in the sphere with the following properties. There is a 
circle corresponding to each face and to each vertex. The vertex circles form a 
packing with two circles touching if and only if the corresponding vertices are 
adjacent. Likewise, the face circles form a packing with circles touching if and 
only if the corresponding faces are adjacent. For each edge, there is a pair of 
touching vertex circles and a pair of touching face circles. These pairs touch in 
the same point, intersecting each other orthogonally. 

This circle pattern is unique up to Mobius transformations. 

If the cellular decomposition is a triangulation, the additional condition is 
automatically fulfilled: Three mutually touching circles have three points of 
contact. There is always a circle through them, and this circle intersects the 
three original circles orthogonally. The above theorem is therefore indeed a 
generalization of Koebe's theorem. 

Furthermore, theorem |l| is equivalent to the following theorem, which "prob- 
ably yields the best and most powerful version of Steinitz' theorem |Zic95| ." Q 



The original Steinitz theorem only asserts the existence of a three dimensional convex 
polyhedron which is combinatorially equivalent to a polytopal cell decomposition of the sphere. 
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Theorem 2. For every polytopal cellular decomposition of the sphere, there is 
a combinatorially equivalent polyhedron with edges tangent to a sphere. This 
polyhedron is unique up to projective transformations which fix the sphere. 

There is a simultaneous realization of the dual polyhedron, such that corre- 
sponding edges of the dual and the original polyhedron touch the sphere in the 
same points and intersect orthogonally. 



The history of theorems [l] and ^ is somewhat unclear. Sachs' survey 



provides valuable information. However, some of the articles in preparation 



which are cited have never appeared in print. Thurston Thv,], who revived the 
interest in circle packings and circle patterns by pointing out the connection to 
'discrete conformal maps' and the theory of 3-manifolds, devised a method to 
construct some kinds of circle patterns and thus prove existence and uni quenes s 
theorems. These methods were adapted by Brightwell and Scheinerman |BS93 |, 



who state and prove theorem |. Schramm jgch92H proves a theorem that is 



much more general than theorem |2[ In it, the sphere is replaced by an arbitrary 
convex body. The proof is not constructive. 

Theorem [j] may be generalized in two directions. First, one may consider 
patterns of circles intersecting not orthogonally, but at arbitrary angles. There 
is a correspondence between such circle patterns in the sphere and polyhedra 
in hyperbolic space: In the Poincare ball model, hyperbolic planes are repre- 
sented by spheres that intersect the unit sphere orthogonally. Hence there is a 
correspondence between circles in the unit sphere and hyperbolic planes. Fur- 
thermore, the intersection angle of two circles equals the intersection angle of 
the corresponding planes. Orientation preserving isometries of hyperbolic space 
correspond to Mobius transformations of the sphere at infinity and vice versa. 
We consider circle patterns that correspond to hyperbolic polyhedra with ver- 
tices in the sphere at infinity, so called 'ideal' hyperbolic polyhedra. These are 



classified by the following theorem, which is due to Rivin Riv96 



Theorem (Rivin). Let £ be a polytopal cellular decomposition of the sphere. 
Suppose an angle 9 e , with < 9 e < ir, is assigned to each edge e. There exists 
an ideal hyperbolic polyhedron which is combinatorially equivalent to S, and 
with exterior dihedral angles 9 e , if and only if the following condition holds: For 
every cocycle 7 of edges, ^2 eEl e > 2n, with equality if and only if 7 is the 
coboundary of a single vertex. 

This ideal hyperbolic polyhedron is unique up to isometry. 

Here, the exterior intersection angle of two faces of a polyhedron is the com- 
plement to 7r of the interior intersection angle, which is the angle inside the 
polyhedron. The exterior intersection angle of two circles is the angle 9 e of fig- 
ure H it is the complement to n of the interior intersection angle 9* . A cocycle 
is a cycle in the Poincare dual cellular decomposition, faces of which correspond 
to vertices of the original decomposition, vertices to faces, and edges to edges. 
The concept dual to the boundary of a face is the coboundary of a vertex. 

Corollary. Let £ be a polytopal cellular decomposition of the sphere, and sup- 
pose every vertex has n > 3 edges. Suppose that every cocycle which is not 
the coboundary of a single vertex is more than n edges long. Then there exists, 
uniquely up to Mobius transformations, a circle pattern on the sphere which is 
combinatorially equivalent to £ and has exterior intersection angles 2ir/n. 
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(a) (b) 

Figure 3: A cellular decomposition (a) and its medial decomposition (b). 

The case n = 4 implies theorem [j], and hence theorem ||. Indeed, suppose 
E is a polytopal cellular decomposition of the sphere. In theorem [l], circles 
correspond to faces and vertices. To apply the corollary, consider the medial 
cellular decomposition of E, whose faces correspond to the faces and vertices 



of E and whose vertices correspond to edges of E. Figure 3(a) shows part of 
a cellular decomposition of the sphere and figure |3(b)| its medial. The dotted 



lines in figure 3(a) represent the edges of the medial decomposition. In |3(b)| , the 



faces of the medial decomposition which correspond to vertices in the original 
decomposition are shaded. To deduce theorem |l| from the corollary, one has to 
show that if we start with the polytopal cellular decomposition E, its medial will 
be polytopal, have 4-valent vertices, and all cocycles which are not coboundaries 
of a single vertex will be longer than 4. In fact, this is easy to see. Suffice it to 
remark that Steinitz' Bedingung des Nichtiibergreifens for E corresponds to the 
property of the medial that cocycles are longer than 4. 

The second direction in which theorem [l] may be generalized is to consider 
not only circle patterns in the sphere but also in other surfaces of constant cur- 
vature. We will consider circle patterns in closed compact surfaces of arbitrary 
genus, as well as circle patterns in surfaces with boundary. 

It is enough to consider only orientable surfaces. To deal with non-orientable 
surfaces, consider first the orientable double cover; the uniqueness of the circle 
patterns in question will imply that the pattern on the cover can be projected 
down to the non-orientable surface. 

By way of a further generalization, one may also consider circle patterns in 
surfaces which may have cone-like singularities in the centers of circles and in 
the intersection points. 



Colin de Verdiere [CdV91] considers Euclidean and hyperbolic circle pack- 



ings, Bragger [Bra92 only the Euclidean case. Rivin [Riv94 and Leibon [LciOl 
treat Euclidean and hyperbolic circle patterns, respectively. They allow cone- 
like singularities in the vertices. 

In this paper, we present new variational principles for circle patterns in Eu- 
clidean surfaces (with curvature 0) and in hyperbolic surfaces (with curvature 
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— 1). The functional are functions of the radii of the circles. Critical points 
correspond to circle patterns with prescribed intersection angles and cone-like 
singularities. Circle patterns in the sphere are treated by projecting them stere- 
ographically to the plane and then using the Euclidean functional. 

We give a complete proof of Rivin's theorem, and we prove an analogous 
theorem for surfaces of higher genus, theorem 14. Both are deduced from the 
more fundamental, but less geometric theorem |f 

The idea to treat circle patterns and cir cle packings using a variational prin- 
ciple is not new. Colin de Verdiere | CdV9l| was the first to present a variational 
principle for circle packings, in which the circles correspond to the vertices of 
a triangulation. He treats both the Euclidean and the hyperbolic case. How- 
ever, only the derivatives of the functionals are given. This suffices to prove the 
existence and uniqueness of circle packings. Bragger [Bra92] remarks that, inves- 



tigating Colin de Verdiere's functionals, one quickly arrives at the Lobachevski 
function. But he does not present these functionals in closed form. Instead, he 
derives a new variational principle. Rivin [Riv94], and Leibon [LciOl present 
variational principles for circle patterns. The circles correspond to the faces of 
a triangulation. Rivin treats the Euclidean, Leibon the hyperbolic case. 

It was unclear how these functionals were related. We show how the function- 
als of Colin de Verdiere, Bragger, and Rivin can be derived from our functionals. 
We believe that Leibon's functional can be identified with a functional we derive 
in section ^|. 

Finally, we would like to add some remarks on the practical side of actu- 
ally constructing circles patterns. The functionals of Colin de Verdiere are not 
presented in closed form; only formulas for their derivatives are given — and in 
a form which seems to have prevented their integration. The functionals of 
Bragger, Rivin, and Leibon are presented in closed form, but their minimiza- 
tion has to be constrained to spaces of 'coherent angle systems', which makes 
their application more difficult. Our functionals are given in closed form and 
there are no constraints. 

Thurston's proof of Koebe's theorem is also constructive. "[It is] based on 
a practical algorithm for actually constructing patterns of circles. The idea of 
the algorithm is to adjust, iteratively, the radii of the circles. A change of any 
single radius affects most str ongly the curvature at that vertex, so this process 

p. 347] Mohar shows that such an algorithm 
. Stephenson's program circlepack constructs 



converges reasonably well ." |Thu 



runs in polynomial time [Moh93 



circle packings by Thurston's method. It was used by Dubejko and Stephenson 
for experiments in discrete analytic function theory [DS95|. There are also 



efforts to apply such circle packing algorithms in medical imaging [HBS + 99, 
MarOl[ . 

Thurston's algorithm can be interpreted as a special method to minimize our 
functionals: Adjusting a single radius so that the neighboring circles fit around 
it is equivalent to minimizing our functionals in one coordinate direction. 



3 Main results 

After some preliminaries, we state the main results, theorem || and theorem |j. 
Theorem ^ is a fundamental (but somewhat technical) existence and uniqueness 
theorem for circle patterns. Section |I] is devoted to the proof of theorem [| From 



G 



Figure 4: A Euclidean kite-shaped quadrilateral, 9 e + 6* = tt. 



it, we deduce (in section [sj) Rivin's theorem and its higher genus analogue, 
theorem |J. 

Suppose a graph is embedded in a closed topological surface such that the 
complement of the graph is the disjoint union of discs. This is called a 2- cell 
embedding of the graph. Such a 2-cell embedding of a graph gives rise to a cellu- 
lar decomposition of the surface. Two decompositions are considered equivalent 
if one is mapped to the other by a homeomorphism of the surface. We only con- 
sider the case of oriented surfaces. A cellular surface is a cellular decomposition 
arising from a 2-cell embedding, considered up to homeomorphism. Associated 
with a cellular surface is a set F of faces, a set E of oriented edges, and a set 
V of vertices. By E we denote the set of unoriented edges. For each unoriented 
edge in E, there is a pair of oppositely oriented edges in E. 

In the case of a surface with boundary, consider a graph embedded in the 
surface, such that the graph meets the boundary only in vertices, and such that 
the complement of the graph and the boundary is a disjoint union of open discs. 
This gives rise to a cellular decomposition. As the edges of the corresponding 
cellular surface we consider only the edges of the embedded graph. This means 
that there are no boundary edges: Every edge has a face on the left and on the 
right. But there are faces whose boundary is not closed, the boundary faces, 
and there are vertices whose coboundary is not closed, the boundary vertices. 

A convenient way to represent a cellular surface E jointly with its dual E* is 
the dual of the medial decomposition, E 9 , which is called a quad-graph, since all 
its faces are quadrilaterals. They correspond to the edges of E. There are two 
kinds of vertices of E 9 , 'white' vertices corresponding to faces of E and 'black' 
vertices corresponding to vertices of E. Q uad-graphs turn out to be crucial for 
defining discrete integrable systems [ BS02[] (in p articular those related to circle 
patterns) and discrete Riemann surfaces Mcr01 |. 

Suppose we are given a positive radius r/ for each face / of E and an angle 
9 e € (0, tt) for each non-oriented edge of E. For each face of E g , construct a 
(Euclidean, spherical, or hyperbolic) kites-shaped quadrilateral as in figure |]. 
Because corresponding sides have the same length, these kites may be glued to 
together to form a cellular decomposition equivalent to E g . One obtains not 
only a cellular decomposition of a topological manifold but a circle pattern on a 
constant curvature surface with cone-like singularities. The circles correspond 
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to faces of E. The vertices of E correspond to points of intersection. 

The cone-like singularities are in the centers of circles and in the points 
corresponding to vertices. The cone angle at an interior vertex v £ Vs is 



e„ = ^0 e , (i) 



where the sum is taken over all oriented edges with initial vertex v. That is, we 
interpret 9 as a function on the oriented edges with 9 e = 9- e . (If an edge e has 
v as initial and terminal vertex, then 9 e appears twice in the sum.) 

The cone angle in the center of a circle corresponding to an interior face / 

is 

$ f =2j2<Pe, (2) 

where the sum is taken over all oriented edges in the oriented boundary of /. 
An elementary geometric calculation yields for the angles 

= 77J js- I 3 ) 



2i rfj - rfk e . 



in the Euclidean case, and 
1 



, tanh f h^-J - tanh (J^) e~ 

Ve 2l \ tanh (frf)- tanh e ; 




in the hyperbolic case. Here, fj and fk are the faces to the left and right of e. 

For boundary faces and vertices, <&/ and 0„, given by the same formulas, are 
the interior angles of the boundary polygon of the quadrilateral decomposition 



(the dotted polygon in figure |2(b)D 



The idea is to prescribe only the intersection angles 9 e and the cone (or 
boundary) angles in the faces, and to solve for the radii. The cone (or 
boundary) angles 0„ in the vertices are determined by the intersection angles 
and equation (0). We will prove the following theorem. 

Theorem 3. Let E be an oriented cellular surface, let 9* G (0,ir) E be a function 
on the non-oriented edges, and let $ G (0, oo) Fs be a function on the faces. 

A Euclidean circle pattern which is combinatorially equivalent to E, which 
has interior intersection angles 9* , and which has cone (or boundary) angles $/ 
in the centers of the circles exists, if and only if the following two conditions 
are satisfied. 



$>(/) = £>*(<0 (5) 
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(ii) If F' is a nonempty subset of the face set F, F' ^ F, and E' is the 
set of all edges incident with any face in F' , then 



£ $(/) < 2f*(e). (6) 

/£F' e£E' 

The Euclidean circle pattern is unique up to similarity. 

A corresponding hyperbolic circle pattern in a surface of constant curvature 
— 1 with cone-like singularities exists, if and only if inequality (j^) holds for all 
nonempty subsets F' C F, including F itself. The hyperbolic circle pattern is 
unique up to isometry. 

This theorem is proved in section || 

Equation (^|) , which is necessary for the existence of a Euclidean circle pat- 
tern, is equivalent to the Gauss-Bonnet formula 

E K f + E K * = 2 ^< 

feF vev 
where x is the Euler characteristic of E, 



Kf 

and 

K„ = 



2ir — $/ for interior faces 
7r — for boundary faces 



2tt — Q v for interior vertices 
7r — Q v for boundary vertices 



Note that the formula \ = \F\ — \E\ + \V\ does not hold for cellular surfaces 
with boundary faces and boundary vertices. To get a topological invariant, one 
must either add boundary edges first, or count the boundary faces and vertices 
only half. 

So far, we did not consider circle patterns in the sphere. The metric ap- 
proach, focusing on the radii of circles, is less useful in this case. The reason 
is this: The only circle preserving transformations of the hyperbolic plane are 
isometries, those of the Euclidean plane are similarities. The Mobius transfor- 
mations, which are the circle preserving transformations of the sphere, distort 
the metric in a more complicated way. 

However, spherical circle patterns without cone-like singularities may be 
treated by projecting them stereographically to the plane. This is described 
in detail in section where, in particular, we give a proof of Rivin's theorem. 

Searching for a generalization of Rivin's theorem to surfaces of higher genus, 
one may conjecture the condition to be that the angle sum for a coboundary is at 
least 2ir, and equal to 2tt only for coboundaries of a single vertex. This is not the 
right condition, however. Firstly, one needs to consider only the coboundaries of 
simply connected domains. Secondly, it is not enough to consider coboundaries. 



For example, if we cut the torus of figure 1(a) along the dual edges represented 
by the heavy dashed lines, it is separated into two pieces; one of these pieces 
is a disc. While it is indeed necessary that the angle sum over the edges we 
cut along is strictly greater than 2n, these edges do not form a boundary. To 
formulate the right condition, we consider cutting the dual cellular surface along 
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some of its edges. Of the pieces into which this separates the surface, only the 
discs have to be considered. The angle sum over their boundary has to be at 
least 2ir. What counts is the boundary of the pieces after cutting, not their 
boundary as regions in the surface. 

Theorem 4. Let S be a cellular surface of genus g > 0, and let T,* be the dual 
cellular surface. Suppose 9 G (0,ir) E is a function on the edges of S, which 
sums to 2tt around each vertex. We identify the edges of £ with the edges of E* 
and consider 9 also as function on the edges o/S*. 

There exists a circle pattern in a surface of constant curvature —1 which is 
combinatorially equivalent to £ and has exterior intersection angles 9, if and 
only if the following condition is satisfied. 

Suppose we cut S* along a subset of its edges, obtaining one or more 
pieces. For any piece that is a disc, the sum of 9 over the edges in 
its boundary has to be at least 2tt, with equality if and only if the piece 
consists of only one face of S* (vertex ofYi). 

Both the constant curvature surface and the circle pattern are determined up to 
isometry if g > 1 and up to similarity if g = 1 . 

The proof is contained in section 



4 The proof of theorem [3| 



In this section, we prove theorem^. Our proof is based on a variational principle 
for circle patterns. The corresponding functional s, on e for the Euclidean and one 



for the hyperbolic case, are presented in section 4T. They are functions of the 
radii. Of course, they also depend on the cellular surface; and the prescribed 
intersection angles and cone angles enter as parameters. A radius function 
which is a critical point yields a circle pattern with the prescribed intersection 
angles and cone angles. The functionals are strictly convex, and this implies the 
uniqueness claims of theorem |J 

Proving the existence claims is a little more involved. One has to estimate the 
functionals from below to show that they become large if some radii tend to zero 



or infinity. In section 4.2, we define 'coherent angle systems', which are functions 
on the directed edges satisfying some linear equations and inequalities. We show 
that the functionals have minima if and only if a coherent angle system exists. 
This reduces the problem to the feasibility problem of linear programming, or 
equivalently, to the feasible flow problem of network theory. 



In section 4.3, applying the feasible flow theorem of network theory, we show 
that the conditions of theorem [| are necessary and sufficient for the existence 
of a coherent angle system. 

4.1 The Functionals. Convexity and the uniqueness of 
circle patterns 

Instead of the radii r, it is more convenient to use transformed variables p, the 
transformation being 



P = 



log r in the Euclidean case 

logtanh(r/2) in the hyperbolic case 
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In the Euclidean case, positive radii correspond to real p. In the hyperbolic case 
case, positive radii correspond to negative p. 

The problem is to find a function p : F — > K (which must be negative in the 
hyperbolic case) satisfying equation (Q) with 

I 1 _ e Ph ~ ?U 'We 

in the Euclidean case and 

1 1 _ e Pfk ~ PS, - ! I _ e Pf k + 9Si - We 
(fie = — log 77. log T7T-. (8) 

* 2* 1 _ e Ph - Ph +W e 2^ S 1 _ ePh + p h +i9 e K > 

in the hyperbolic case. Here, fj and are the faces to the left and right of e, 
respectively. (Equations (Q) and (||) are equivalent to equations @ and (jl]).) 

We turn this into a variational problem. The corresponding functionals are 
defined in terms of the imaginary part of the dilogarithm function Li2. In 
the appendix, all requisite facts about this function and the related Clausen's 
integral CI2 are summarized. 

Definition. The Euclidean circle pattern functional is 

Seuc(p) = £ (lmLi 2 [ePh ~ Pf> + + ImLi 2 (ePf* - Ph + 

-K{pti ,+Ph)) +J2 (p fpf- 

of 

The first sum is taken over all edges e, and fj and fk are the faces on either side 
of e. The summands are symmetric in fj and fk, so it does not matter which 
face is considered as fj and which as fk- The second sum is taken over all faces 
/■ 

The hyperbolic circle pattern functional is, using the same notation, 
Shy P (p) = £ (im Li 2 (ePh ~ Ph + # e) + Im Li2 ( e P/, - Ph + Me) 

+ ImLi 2 (e^ + ^ + + ImLi 2 (e - ^ ~ ^ + 

(10) 

Proposition 1. The critical points of the functionals Se U c and Sh yp correspond 
circle patterns with intersection angles 6 and cone angles <E> in the centers of 
circles. In particular, if p is a critical point of Sh yp then all p / are negative. 

Proof. The functionals are of the form 
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where F- is an even function. The partial derivative with respect to a face pf j 
is 

H = E (*"-(Pf>-Ph)+no>f i +Pfu))+G'o>f i ), 

where the sum is taken over all oriented edges in the boundary of fj. Since, 
from equation (^q) of the appendix, 

^-ImLi 2 ( e ^ + ^) = ilog 1 -' 



dx ' 2i i_ e xl» 

one finds that 

^- = -2(E^e)+%, (12) 

where the sum is taken over all directed edges in the boundary of /. 

It is left to show that if p is a critical point of Shyp, then all pf are negative. 
The function 

1 1 - e x ~ ^ e 
fW^-log——^- 

is strictly increasing for < 9 < it. Hence, if, in equation (||), pf. is nonnega- 
tive, then <£e is non-positive. Therefore, if some pf is nonnegative, then, from 
equation (|T4), 

^>$(/)>0, if Pf >0. (13) 
dp/ 

□ 

Rescaling all radii by a factor h corresponds to the transformation p i— > 
/? + /i If, where If is the function that is constantly 1 on F. Note that 

S E uc(P + hl F ) = SEuc(p) + h(j2®f- 2 J2 9 e 

fEF e£E 

The functional can have a critical point only if the coefficient of h vanishes. In 
this case, the functional is scale-invariant, and one may restrict minimization to 
the subspace 

U = { P eR F \Y / Pf = 0}- (14) 
feF 



Proposition 2. The hyperbolic functional is strictly convex. The Euclidean 
functional is strictly convex on the subspace U. This implies, in the hyperbolic 
case, the uniqueness of circle patterns with prescribed intersection angles and 
cone angles, and, in the Euclidean case, uniqueness up to scale. 
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Proof. Taking second derivatives on both sides of equation (|Tlh, one obtains 



S " = J2 (F- (P fj ~ Pfk ) id Pfi - d Pfk f + F» ( Pfj + Pfk ) {d Ph + d PSk f 



+ Y / G"(p f )dp}. 

of 



Thus, by a straightforward calculation, 



SLc = £ ooah(p A - P/ ;)-ooBfl.^-^ : 



and 



(d P/ . +dp fk ) 



cosh( / j / . +p /fc ) -cos( 

These second derivatives are positive definite on U and M F , respectively. □ 

Remark. Formulas (^) and ( JlO| ) for 5e uc and SJ^p involve the imaginary part of 
the dilogarithm of a complex argument. The formulas below involve only real 
functions. Since, moreover, Clausen's integral is periodic, these formulas are 
well suited for the practical evaluation of the functionals. 

For p G K F , let p and s be the functions on oriented edges defined by 

p e = 2 arctan [tan ( ) tanh ( Pfk - Pfl ) ) (15) 



2 



and 



s e = 2 arctan ^tan tanh y^- JJ > ( 16 ) 

where /j and are the faces to the left and right of e. Note that s e = s_ e (so 
that s may be considered as a function on non-oriented edges), but p e = — p_ e . 
The following formulas hold: 



&„ c (p) - ^ (p e ( PA - p,.) + Ch(9* e + Pe ) + Ch(0* e -p e ) - Cl 2 (2^) 

~ KiPfj + Pf k j) + , 
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and 

Shyp(p) = E (PA "/>/*)+ C1 2 ( e + Pe) + Cl 2 (fl e * - Pe ) 

+ Se ( P/j + Pfk ) + c\ 2 (o* e + s e ) + ci 2 (e; - s e ) 

-2C\ 2 {26$) +Y,*fPf- 

°f 

The first sum in each equation is taken over non-oriented edges, or rather, 
over one oriented representative e for each non-oriented edge. Since both p e 
and (fk — fj) change sign when — e is chosen instead of e, the summands are 
independent of which choice is made. 

These formulas can be derived using equation ( fio|) of the appendix and some 
simple properties of Clausen's integral, which are also found there. 

The hyperbolic functional will reappear in this form in equation ([35]) of 
section p. 



4.2 Coherent angle systems and the existence of circle 
patterns 



In section 4.1, the uniqueness of a circle pattern was deduced from the convexity 
of the functionals. The rest of section ^ is devoted to the existence part of 
theorem |^. To establish that the functionals have a minimum, we show that 
they tend to infinity as p tends to infinity in R or U, respectively. 

The main difficulty in estimating the functionals from below is to somehow 
merge the sum over edges with the sum over faces in equations (^) and (|To|). 
This is achieved with the help of a so called 'coherent angle system'. Therefore, 
in this section we prove that the functionals have minima if and only if coherent 
angle systems exist. In section 4.3, we show that conditions of theorem [jj are 



necessary and sufficient for the existence of a coherent angle systems. 

Definition. A coherent angle system is a function ip £ M. E on the oriented edges 
satisfying, in the Euclidean case, 

1. for all oriented edges e G E, 

tp e > and tp e + tp_ e = 9* e 

2. for all oriented faces / £ F, 

where the sum is taken over all oriented edges in the oriented boundary 
of/. 

In the hyperbolic case, condition [j] has to be changed to 
|l|'. for all oriented edges e € E, 

ip e > and tp e + tp- e < 9*. 
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Proposition 3. The functional Se U c (Shyp) has a critical point, if and only if 
a Euclidean (hyperbolic) coherent angle system exists. 

Proof. If the functional Se U c {Shyp) has a critical point p, then equation (Q) 
(equation (g)) defines a coherent angle system. It is left to show that, conversely, 
the existence of a coherent angle system implies the existence of a critical point. 

Consider first the Euclidean case. Suppose a Euclidean coherent angle sys- 
tem <p exists. One easily deduces equation (|^). Hence, the functional Se U c is 
scale invariant. We will show that Se U c(p) oo if p ^ oo in the subspace U 
defined in equation More precisely, we will show that for p € U, 

Smc(p) > 2mhxip e max (17) 

eeE f^ F 

The functional Se U c must therefore attain a minimum, which is a critical point. 
For x G K and < 6 < ir, 

ImLi 2 (e- T + id ) + ImU 2 {e- x + id ) > (tt - 9) \x\, 

and hence 

Seuc(p) > -2 9 e min iPh ' Pf,) + I] */P/. 

Throughout this section, '5Z egS ' ^ s mean t to indicate that the sum is taken 
over one representative from each pair of oppositely oriented edges; it does not 
matter which. The faces on either side of e are fj and fk- 

Now we use the coherent angle system to merge the two sums. Because 

®fPf = 2 J2( ( Pe Pfj + V>-e Ph), 
f£F eeE 



one obtains 



Seuc(p) > 2 ^ min (<Pe, <P-e) \pf k ~ Ph I 



eeE 

Since we assume the cellular surface to be connected, we get 

Seuc(p) > 2min<p e (max/)/ - minp/), 
eeE K feF rj feF tJ > 

and from this the estimate (|l7j). 

The hyperbolic case is similar. One shows that, if all pf < 0, 

Shyp{p) > 2min|<p e + tp-e - 0*1 max \pr\. 
yF eeE 1 1 feF 1 1 



Since it was shown in section 4J that dSh yp /dpf > if pf > 0, this implies the 



existence of a minimum. □ 

4.3 Conditions for the existence of coherent angle systems 

With this section we complete the proof of theorem |j. All that is left to show 
is the following proposition. 
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Figure 5: The network (iV, X). Only a few of the branches and capacity intervals 
are shown. 



Proposition 4. A coherent angle system exists if and only if the conditions of 
theorem Q hold. 

Proof. It is easy to see that these conditions arc necessary. To prove that they 
are sufficient, we apply the feasible flow theorem of network theory. Let (N, X) 
be a network (i.e. a directed graph), where N is the set of nodes and X is the 
set of branches. For any subset N' C N let ex(N') be the set of branches 
having their initial node in N' but not their terminal node. Let in(N') be the 
set of branches having their terminal node in N' but not their initial node. 
Assume that there is a lower capacity bound a x and an upper capacity bound 
b x associated with each branch x, with — oo < a x < b x < oo. 

Definition. A feasible flow is a function ip G M. x , such that Kirchoff's current 
law is satisfied, i.e., for each n E N, 

xG ex({n}) x£in({n}) 

and a x < ip x < b x for all branches x. 

Feasible Flow Theorem. A feasible flow exists, if and only if for every nonempty 
subset N' (Z N of nodes with N' ^ N, 

Y bx - Y ax - 

x£ex(N') xein(N') 



A proof is given by Ford and Fulkerson [FF62, Ch. II, §3]. (Ford and Fulk- 
erson assume the capacity bounds to be non- negative, but this is not essential.) 

To prove the proposition in the Euclidean case, consider the following net- 
work. (See figure ||). The nodes are all faces and non-oriented edges of the 
cellular surface, and one further node that we denote by M: N = F U E U {Kl}. 
There is a branch in X going from M to each face / G F with capacity interval 
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[i<i>/,oo). From each face / there is a branch in X going to the non-oriented 
edges of the boundary of / with capacity interval [e, oo), where e > will be 
determined later. Finally there is a branch in X going from each non-oriented 
edge e e E to S3 with capacity (— oo, 0*]. 

A feasible flow in the network yields a coherent angle system. Indeed, Kir- 
choff's current law at S3 implies that the flow into each face / is h$f and the 
flow out of each edge e is #* . It follows that the flow in the branches from F to 
E constitutes a coherent angle system. 

Assume the conditions of theorem H] are fulfilled. We need to show that the 
condition of the feasible flow theorem is satisfied. Suppose N' is a nonempty 
proper subset of N. Let F' = N' n F and E' = N'C\E. 

Consider first the case that S3 6 N', which is the easy one. Since N' is a 
proper subset of N there is a face / S F or an edge e € E which is not in N'. In 
the first case there is a branch out of N' with infinite upper capacity bound. In 
the second case there is a branch into N' with negative infinite lower capacity 
bound. Either way, the condition of the feasible flow theorem is trivially fulfilled. 

Now consider the case that S3 $ N'. We may assume that for each face 
feF', the edges of the boundary of E are in E'. Otherwise there are branches 
out of N' with infinite upper capacity bound. For subsets A,BcN denote by 
A — > B the set of branches in X having initial node in A and terminal node in 
B. Then the condition of the feasible flow theorem is 

£ U f + e\F\F'^E'\< Y^K- 

feF' ee£' 

It is fulfilled if we choose 




where the minimum is taken over all proper nonempty subsets F' of F and 
E'(F') is the set of all non-oriented edges incident with a face in F 1 . The 
minimum is greater than zero because of condition (ii) of theorem ^. 

In the hyperbolic case, the proof is only a little bit more complicated. In the 
network, the flow in the branches going from S3 to a face / must constrained to 
be exactly |<&(/); and the capacity interval of branches going from an edge e 
to S3 has to be changed to (— oo, 6*(e) — e]. □ 

5 Circle patterns in the sphere, Rivin's theo- 
rem, and its higher genus analogue 

In section |5.l| we adapt theorem || to circ le patter ns in the sphere. Then we use 



it to prove Rivin's theorem in section |5.2| . Section |5.3| is devoted to the proof of 
theorem |[ the higher genus analogue. 

5.1 Circle patterns in the sphere 

A circle pattern in the sphere may be projected stereographically to the plane, 
choosing some vertex, Voo, as the center of projection. For example, figure 
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(a) (b) 

Figure 6: (a) The regular cubical pattern after stereographic projection, (b) 
How to calculate the $ / for the new boundary circles. 



3(a) shows the circle pattern combinatorially equivalent to the cube and with 
intersection angles 7r/3 after stereographic projection. 

One obtains a circle pattern in the plane in which some circles (those cor- 
responding to faces incident with Voo) have degenerated to straight lines. Since 
stereographic projection is conformal, the intersection angles are the same. Fur- 
thermore, a Mobius-equivalent circle pattern in the sphere leads to a planar 
pattern which is similar, provided the same vertex is chosen as the center of 
projection. 

The idea is to construct a circle pattern in the sphere by constructing the 
corresponding planar pattern using the Euclidean functional and then projecting 
it to the sphere. 

To construct the planar pattern, we first remove from the cellular surface the 
faces incident with and all edges incident with them. Some of the remaining 
faces will then be boundary faces. For such a new boundary face / we set 

^ f = 2n-J2^l (19) 



where the sum is taken over all removed edges incident with /; see figure |6(b) 
For all interior faces /, we set 

$ / = 27r. (20) 

Assuming the conditions of theorem ^| are satisfied, we then construct the 
corresponding planar pattern, add the lines corresponding to the removed faces, 
and project to the sphere. 

In this way, we deduce the following proposition from theorem ^. 

Proposition 5. Let £ be a closed cellular surface of genus 0, and let an angle 
be assigned to each non-oriented edge by 9 E (0,ir) E . Suppose that for each 
vertex v E V , the sum of 9 over the edges around v is 2tt, so that there is no 
curvature in the vertices. Let 9* = n — 9. 

Suppose E V, and let Foq be the set of all faces incident with v^. Suppose 
the following conditions holds. 
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(i) The 1-skeleton of the Poincare dual £* remains connected if the faces 
in Frx, and all edges incident with them are removed. 

(ii) 

2n(\F\-\F 00 \) = J2^*(e), (21) 

where the sum is taken over all edges incident with a face in F\F OOJ i. e. 
all edges not incident with Voo . 

(Hi) If F' is a nonempty subset of F \ Foo and F' ^ F \ F^, and E' is 
the set of all edges incident with a face in F' , then 

2w\F'\ < J2 26l *( e )- ( 22 ) 

Then there exists a circle pattern in the sphere which is combinatorially equiva- 
lent to X and which has intersection angles 9* . This circle pattern is unique up 
to Mobius transformations. 

Proof. By condition (i), after we remove from £ all faces in Foo and all edges 
incident with them, we are left either with a single face, or with a connected 
cellular surface Eo with boundary faces and boundary vertices. The first case 
is elementary. In the second case, we want to apply theorem ^|. 

For boundary and interior faces / of So define $/ by equations (|l^) and 
(pp[), respectively. We have to show that for boundary faces $/ > 0. Suppose 
/ is a boundary face of So- Let F' — F \ (F^ U {/}). Subtracting ( pi] ) and 
( [22|) , one obtains 2tt > ^2 2#* , where the sum is taken over all removed edges 
incident with /. Hence $/ > 0. 

We apply theorem || to the cellular surface So with intersection angles given 
by 9 restricted to the edges of E and with $ as defined above. Conditions (ii) 
and (iii) of the proposition imply conditions (i) and (ii) of theorem |[ Indeed, 
the face set of £0 is F \ F^ and 

$/ = 2 7 r(|F|-|F 00 |)-^2C 

f&F\Foo 

where the sum on the right hand side is taken over all edges of S that are not 
edges of £0 but are incident with a face in F \ F^ . Hence, condition (ii) of the 
proposition implies 

feF\F x 

where the sum is taken over all edges of £oj i- e - condition (i) of theorem 

Similarly, condition (iii) of the proposition implies condition (ii) of theo- 
rem [}[ Therefore, the corresponding planar circle pattern exists, and by the 
construction described above, we can construct the spherical circle pattern. □ 
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Figure 7: The circles and quadrilaterals corresponding to a simple closed copath 
after stereographic projection. 

5.2 A proof of Rivin's theorem 

We now give a proof of Rivin's theorem. First, we show that the condition of 
Rivin's theorem is necessary. If 7 is the coboundary of a single vertex, then 
X)e£7 = 27T simply means that the corresponding circle pattern has no cur- 
vature at that vertex. Suppose a circle pattern on the sphere is given and 7 is 
a closed copath that does not bound a single vertex. Without loss of generality, 
we may assume that 7 is simple, i. e. that it meets a face at most once. 

Suppose that 7 does not meet all the faces of the cellular decomposition; we 
will deal with the other case afterwards. Then there is a point in the sphere 
which lies outside of all the circles belonging to 7. (In fact, each circle of 
a circle pattern contains points which are not contained in any of the other 
circles.) Project the circle pattern stereographically to the plane, using such a 
point as the center of projection. One obtains a sequence of intersecting circles 
as in figure [?|. Now 

E(*-«i-&-7j) = 27r (23) 

and 

8 J= n- a : - pi (24) 

This implies ^2 ®j — 27T, with equality if and only if all 7^ are zero. This is the 
case only if the copath bounds a single vertex. 

Now suppose that the copath 7 does meet all the faces of the cellular decom- 
position. Consider the graph (V, E\j), consisting of all vertices of the cellular 
decomposition and all edges that do not appear in 7. This graph is the disjoint 
union of two trees. Therefore, \V\ — \E\j\ — 2. Since 8 sums to 2ir around 
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each vertex, 

2n\V\ = J2 2de - ( 25 ) 

eG-E 

With Q e < it one obtains, 

20 e = 2n\V\ -J2 29 e> 2 <\ V \ ~ \ E \ 

7 E\ 7 

and hence the desired inequality. Hence the condition of Rivin's theorem is 
necessary for the existence of a circle pattern. 

Now, using proposition ^, we show that the condition of Rivin's theorem is 
sufficient for the existence of a circle pattern. Thus we assume that X)eG7 ^ e — 
27r for every copath 7, with equality if and only if 7 is the the coboundary of a 
single vertex. We need to deduce the conditions of proposition |5[ Condition (i) 
follows from the assumption that X is polytopal. 

To deduce conditions (ii) and (hi), suppose that F' is a subset of F\F QO , and 
E' is the set of all edges of £ incident with a face in F' . Consider F" = F\F', 
the complement of F', and E" = E \ E' , the complement of E' . For each 
edge in E" , the faces on either side are in F". Now consider the dual cellular 
decomposition £*, and in it, the graph V = (F",E") with vertex set F" and 
edge set E" . As for any graph, we have 

\E"\ - \F"\ =c-n, 

where n is the number of connected components of T and c is the dimension of 
the cycle space. Since the graph is embedded in £*, a cellular decomposition of 
the sphere, we have 

c = r — 1, 

where r is the number of regions into which T separates £*. Since E" contains 
the edges of £ incident with Voo, or, dually, the edges of £* in the boundary of 
Voo , the number of regions is at least two. 

The boundary of each region is a copath in E, and hence the sum of 9 over 
each boundary is at least 2ir. Sum over all boundaries to obtain 

2?rr < ^2 20 e - 

e£E" 

Indeed, each edge in E" appears in at most 2 boundaries. Equality holds if and 
only if every edge is contained in a boundary and each region contains only one 
vertex or all vertices but one. This means equality holds if and only if E" = E or 
E" = d*v for some vertex v; equivalently, if and only if F 1 = or F' = F\ i 7 ^. 
Thus, we have shown that 

2n(\E"\ - \F"\) < 2e ? - 27r (" + !). ( 26 ) 

eS-E" 

with equality if and only if F' = or F' = F \ Fqo . 
Equation (|||) and Euler's formula, 

\F\-\E\ + \V\ = 2, 
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imply 

2ic(\F\ - \E\) = -^26» e + 47T. (27) 

eeE 

Add equation ( pTj ) and inequality ( p6| ) to obtain 

2tt(|F'| - < - 2d e - 27r(n - 1), 

e&E' 

or, because 9* = tt — 9, 

2w\F'\ < Y 26>*+27r(l-n). 

eGE' 

Equality holds if and only if F' = or F' = F \ . Since n > 1, and n = 1 if 
F' = F \ Foe , we have deduced conditions (ii) and (iii) of proposition |[ This 
completes the proof of Rivin's theorem. 



5.3 Higher genus. Proof of theorem || 

In this section, proceeding in a similar way as in the last section, we derive 
theorem |] from theorem |^. 

In one way, the higher genus case is simpler than the spherical case. We 
do not have to transfer the problem to the plane by a stereographic projection. 
We may use theorem [| directly. For this reason, we do not have to assume the 
cellular surface to be polytopal. 

In another way this case is more complicated. When a graph is embedded 
in a surface of genus g > 1, the dimension c of its cycle space and the number 
r of regions into which the graph separates the surface are not related by the 
simple formula c = r — 1. Instead, the following proposition holds. 

Proposition 6. Suppose a non-empty graph T — [V, E) is embedded in a closed 
surface of genus g and separates the surface into r regions. Then 

r 

r-\E\ + \V\ = 2-2g + J2hj, (28) 

i=i 

where hj is the dimension of the first ^-homology group of the j th region. 

(The graph T is non-empty if V 0.) If all regions are simply connected, 
then all hj = 0, and equation ( p8| ) is Euler's equation for a surface of genus g. A 
multiply connected region can be turned into a simply connected one by adding 
hj edges which do not separate the region. A formal proof of proposition O is 



given by Giblin [Gib77, ch. 9]. It is an application of the long exact sequence of 
relative homology and Lefschetz duality. 

First, we show that the condition of theorem^ is necessary for the existence 
of a circle pattern. For g — 1, the flat case, this follows from equations (^3|) and 
( p4[ ) as in the preceding section. (See figure |?].) For g > 1, one has instead 

E(tt - aj - fa - jj) = 2tt + A 
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and 



d 3 = 7r -A J ~a J -f3 J , 

where A is the area enclosed by 7 and Aj is the area of the triangle with angles 
ctj, (3j and 6j. Hence, ^ Oj = 2tt + X^7j + ^ — E Ah an d therefore ^ 6j > 2ix 
with equality if and only if 7 is the coboundary of a single vertex. 

In the following, we show that the condition of theorem ^ is also sufficient 
for the existence of a circle pattern. Suppose the conditions of the theorem [| 
are satisfied. We deduce the conditions of theorem |^. Thus, let F' be a non- 
empty subset of F and let E' be the set of edges incident with any face in F 1 . 
Let F" = F \ F', E" = E\ E', and consider the graph T = (F",E") in the 
dual cellular surface £*. This graph is non-empty, if F 1 ^ F. In that case, 
proposition ^ implies 

\E"\-\F"\=2g-2 + r-J2h j . (29) 

3=1 

Suppose we cut the dual surface S* along the (dual) edges in E" . This would 
separate the surface into pieces. By the boundary of one such piece, we mean the 
boundary after having cut the surface, not the boundary of the corresponding 
region in the surface. For example, if and edge e £ E" has the same region on 
either side, the boundary of that piece contains e twice. Thus, every edge of 
E" appears exactly twice in the boundaries of the pieces into which the surface 
separates. Hence, if jj is the boundary of the j th region, 

j = l ^e£7j ' eG-B" 

By the condition of theorem ^, for a simply connected piece, i. e. for one with 
hj = 0, the sum of 9 over its boundary is at least 2tt, and equal to 2tt only if 
the piece consists of a single vertex. Hence, we have 

r 

2n(r-^ f h j ) < ^ 26 e . 

j=l e£fi" 

Equality could hold only if each region consisted of a single vertex. But this 
would imply E" = E, or E 1 = 0, contradicting the assumption that F' was 
non-empty. If F" is non-empty, equation ( p9| ) implies 

2n(\E"\ - \F"\) < 2n(2g - 2) + ^ 26 e . (30) 

e£E" 

If F" is empty (and hence also E"), this inequality is trivially true if g > 1, and 
equality holds if g = 1. 

Because there is no curvature in the vertices, equation (|25| ) holds here as 
well. With Euler's formula, 

\F\-\E\ + \V\ = 2-2g, 
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this implies 

27r(|F|-|£:|) = 27r(2-2.g)-^2# e . 

e£-E 

Add this equation to inequality ( |30| ) to obtain 

2tt(|F'| - \E'\) <-J2 2^e, 

with equality if and only if F' = F and g = 1. Since 9* = n — 9, this inequality 
is equivalent to 

2tt|F'| < ]T 20 e *. 
This completes the proof of theorem |j. 

6 Other variational principles 

In this section, we derive the variational principles of Colin de Verdiere, Bragger, 
and Rivin from ours. For circle packings, Colin de Verdiere's functionals are 
obtained from our functionals Se U c and Sh yv for orthogonally intersecting circles 
by minimizing with respect to the radii of those circles which do not appear 
in the pattern. The derivation of Bragger's and Rivin's functionals involves a 
Legendre transformation of Se U c- We believe that Leibon's functional can be 
derived in the same way from Sh yp by the Legendre transformation which we 
present. 

6.1 Colin de Verdiere's functionals 

Colin de Verdiere considers circle packings in which the circles correspond to 
the vertices of a triangulation. He considers the 1-form 

LL> = adu + Pdv + j dw 

on the space of Euclidean triangles, where 

u = logic, t> = logy, and w — logz 

and x, y, z and a, /?, 7 are as shown in figure ||. 

It turns out that dw — 0, hence one may integrate. Define the function 

f ao ,/3 ono on R 3 by 

p(u,v,w) 

ia ,/3 ,7o( u > u ' w ) = ( a o ~ ") du + (/3 - f$) dv + (70 - 7) dw. 

The initial point of the integration does not matter. 

Suppose we are given a triangulation and a coherent angle system for it. 
Here, a coherent angle system is a positive function on the set of angles of the 
triangles, such that the sum in each triangle is it, and the sum around each vertex 
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Figure 8: A Euclidean triangle 



is 2it. For a function p on the vertices of the triangulation, Colin de Verdiere's 
functional for Euclidean circle packings is 



where the sum is taken over all triangles, u, v, and w are the values of p on the 
vertices of each triangle, and ao, Po and 70 are the corresponding angles of the 
coherent angle system. 

The critical points of this functional correspond to the logarithmic radii of 
a circle packing. 

The hyperbolic case is treated in the same way, except that now 

u = logtanh(x/2), v — logtanh(y/2), and w ~ logtanh(z/2), 

and in the definition of a coherent angle system, it is required that the sum of 
the angles in a triangle is less than tt. 

To treat circle packings with our functionals, we consider circle patterns with 
orthogonally intersecting circles on the medial decomposition of the triangula- 



Suppose So is a triangulation. Let £ be the medial decomposition. Set 
9 e = 7r/2 for all edges e. We will consider only the case of closed surfaces, so let 
$/ = 2ir for all faces /. The faces of £ are of two types: those that correspond 
to faces of So and those that correspond to vertices of Eo- Let F = F% U F2, 
where F% contains the faces of the first type and F2 contains the faces of the 
second type. 

Consider S Eu c as function on R Fl x R F2 and define S 2 : R F2 -> K, 



Hence, S^/^) = Seuc{pi, P2), where p\ € M Fl is determined as follows. 
Suppose / 6 F\ and f a , fb, fc S F2 are the neighboring faces of /. Construct 
the Euclidean triangle whose sides are 




tion (see figure ||) . 



S% [pi) 



min S Eu c(pi,P2)- 



e P2f a + e P2f b f e P2f b + e P2f c ; and e P2 /c + e P2f a 



Let pi(f) be the logarithmic radius of the inscribed circle. 
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For p2 £ M , and p\ the corresponding point in 

dSsuc i \ „ OSeuc , \ 9S2 , . 

-3 (Pl,P2)=0, — (P1./92) = -3— (P2j 

dpi dp 2 Cp 2 



Consider S'cdy as function on M 2 . It is not hard to see that dS*2 = dScdv- 
This implies that Colin de Verdiere's Euclidean functional is, up to an additive 
constant, equal to 5 2 - 

Colin de Verdiere's functional for hyperbolic circle packings can be derived 
from Shyp in the same way. 

6.2 A Legendre transformation of the Euclidean functional 

Consider a function 5* : M. F — > M of the form 

S( P )= ]T F e ( Pfk - Pfj ) + J2Gf(Pf^ 
fi°l°fk 0/3 

where the first sum is taken over all oriented edges e and the second sum is 
taken over all faces /. The functional Se U c is of this form with 

F e (x) =ImLi 2 (e x + i9 e) 

and 

G f (x) = ($ f 

where the sum is taken over all edges in the boundary of /. 
Define the function S L : R g x R F -> K by 

S L (v lP )= Y, F M+H G fAPh)- 
fi°l°fk ofj 

The 'L' in the superscript stands for 'Lagrange'. We do not assume that v is 
antisymmetric, i. e. that v e — u_ e . Clearly, the critical points of S correspond 
to the critical points of S L under the constraints 

where fj and fk are the faces to the left and to the right of e, respectively. 

Definition. Suppose F is a smooth and strictly convex real function on some 
open interval. Then 

y = F'(x) (31) 
defines a smooth coordinate transformation. The Legendre transform of F is 

F(y) =xy- F(x), 
where x is related to y by (J3l|). 
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Note that 

p t dF dx ^ ^ dF dx 

dy dy dx dy 

It follows that F is the Legendre transform of F. 

Now let F e be the Legendre transform of F e and define the function S H by 

S H (<p,p)= ]T (MPh ~ Pfi) ~ Fe(<Pe)) +^G( Pfj ). 

The H in the superscript stands for 'Hamilton'. If tp and v are related by 
(fie = F' e (v e ), then S H (ip,p) = S L (v,p). Remarkably the critical points of S 
correspond to critical points of S H , where p and <p may be varied independently: 

Proposition 7. // (}p,p) is a critical point of S H , then p is a critical point of 
S. 

Conversely, suppose p is a critical point of S. Let ip e — F' e {pf k — p/j), where 
fj and fk are the faces to the left and right of e. Then (if,p) is a critical point 
ofS H . 

Proof. We have dS H /d<p e = pf k — pj k — F[,((p e ). So all partial derivatives 
with respect to the variables (p e vanish, if and only if ip and p are related 
by Pfk ~ Pfk = Kife), or, equivalently, tp e = F' e (p h - p h ). In this case 
S H (<p,p) = S(p). □ 

Now consider in particular the Euclidean functional Se U c- The variables p 
and (p are related by equation (Q). The inverse relation is 

P/ fc -^=log sin( ^ + , e) . (32) 



For the functional in Hamiltonian form, identity (39) of the appendix implies 



tfep)= E (ve(p fk -Pf j ) + l{Ch(2Ve)-Ch(2<Pe+26 e )+Ch(26 e )) 



2 

fj°^°fk 



of eedf 



Pf- 



Note that S^ uc depends linearly on p. Collecting the coefficients of each 
Pf, one finds that if is a coherent angle system, then S^ uc (<p,p) does not 
depend on p at all. Therefore, this leads to a new variational principle, where 
the functional depends only on the angles, and the variation is constrained to 
the space of coherent angle systems. 

Proposition 8. If (p is a coherent angle system, then S^ uc (ip,p) — S^ uc (ip), 
where 

SE U M = E ( C1 2( 2 ^) + \ C M20e)) (33) 
and the sum is taken over all oriented edges e. 
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If (p is a critical point of Sjjf nr under variations in the space of coherent angle 
systems, then the equations nSM ) form a consistent system of linear equations 
for the pf. They define a function p S M. F up to an additive constant. This p 
is a critical point of Seuc- 

Conversely, any critical point of Seuc yields, via equation (Qj, a coherent 
angle system which is a critical point of S§ uc under variations in the space of 
coherent angle systems. 



Proof. To derive equation (J33|), note that, if ip is a coherent angle system, then 
<p e + 6 e = 7T - tp_ e . Hence C\ 2 {2<p e + 20 e ) = - Cl 2 (2</>_ e ). 

Now suppose ip is a critical point of S^ uc under variations in the space of 
coherent angle systems. The partial derivatives of S^ uc are 

-^ = -21og(2sin^ e ). 

0<p e 

The tangent space to the space of coherent angle systems is spanned by vectors 
of the form 



d d 



where 7 is some cocycle. (If we add to (p e we have to subtract the same amount 
from p- e . But the sum of ip around each face has to remain constant, so we 
step from one face to the next, adding and subtracting along some cocycle.) 

Since tp is a critical point under variations in the space of coherent angle 
systems, 

v->, sinipe 
P Og sim> e + e ) = ° 



for all cocycles 7. Hence, equations (|32|) are compatible. It is now easy to 
complete the proof. □ 

6.3 Bragger's functional and Rivin's functional 

Like Colin de Verdiere, Br agger | Bra92f considers circle packings in which the 



circles correspond to the vertices of a triangulation. A coherent angle system in 



the sense of section 3T is clearly equivalent to a coherent angle system in the 
sense of this article. Thus, Bragger's functional is seen to be equal to S^ uc /2 up 
to an additive constant. 



Rivin [Riv94] considers Euclidean circle patterns with arbitrary prescribed 
intersection angles. The pattern of intersection is determined by a triangulation. 
Circles correspond to faces of the triangulation. However, since Rivin allows 
intersection angles 8 = 0, adjacent triangles may in effect belong to the same 
circle. Thus, cellular decompositions with non-triangular cases may be treated 
by first dissecting all faces into triangles. Rivin treats cone-like singularities 
in the vertices but not in the centers of the circles. His functional is up to an 
additive constant equal to S§ uc /2. 
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6.4 A Legendre transformation of the hyperbolic Func- 
tional 

Now, consider a function S : M. F — > K of the form 

S (P) = E ( F e(Ph Ph) + F e + ( Pf] + P/J) + £ 

Here, the first sum is taken over all norcoriented edges e, and fj and /fc are the 
faces on either side of e. We assume that F~ is an even function so that is does 
not matter which face is which. The second sum is taken over all faces /. The 
functional Sh yp is of this form with 

F-{x) = F+(x) = ImLi 2 (e x + + ImLi 2 (e~ x + 

and 

Gf(x) = $ f x. 

The space M. E of functions on non-oriented edges splits into the space of 
antisymmetric and symmetric functions, 

= Alt(-E) © Sym(.E), 

where 

A\t(E) = {veR S \ w_ e = -v e } 

and 

Sym(Jj) = {wef fi w- e = w e }. 
Define the function S L : Alt(E) x Sym(£;) x R F -> R by 

Suppose F~ and F + are smooth and strictly convex. The coordinates con- 
jugate to v e and w e are p e = F~'(v e ) and s e — F + '(w e ). Clearly, p_ e = — p e 
and s_ e = s e . Let F~ and F + be the Legendre transforms of F~ and F + . For 
p and s in suitable subsets of Alt(E) and Sym(_E), define 

S H (p,S, P )= {PeiPh- Pf ] )-Fe{p)+Se{p fk +Pf ] )~F+{ S )) 

+£%(/>/)■ 

If p and s are related to p by 

Pe = F?(Pf>-Pf,)> ^ = F+'(p fk +p fi ) ) (34) 

where /j and /& are the face to the left and right of e, then S H (p, s, p) = S(p). 
Again, the critical points of S correspond to critical points of S H , where p, s, 
and p vary independently The following proposition is proved in the same way 
as proposition |7| 
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Proposition 9. If (p, s, p) is a critical point of S , subject only to the con- 
straints p_ e = — p e and S- e — s e , then p is a critical point of S. 

Conversely, suppose p is a critical point of S. Define p and s by equations 
Then (p,s,p) is a critical point of S H . 



Now consider the functional Sh yp . The variables p and s are related to p by 
equations (fl5| ) and ([l6|). The Hamiltonian form of the functional follows from 
equation (HQ) of the appendix: 



s h H yp ( P , s, P )= J2 [Pe (pa -pu) + ci 2 (0; + Pe ) + ci 2 (e; - Pe 

+ **{PU +Ph)+ Cl 2 (0: + s e ) + Cl 2 (8* e - s e ) 



2Ch(26* e ))+J2*fPf- (35) 

of 



(See equation ( [40]) of the appendix.) 
Introduce the variables 



<Pe = 7}(Pe - Se) 

on oriented edges instead of (p, s): 

S h H yp(^P,P) = S" yp (p,S,p). 

This functional is linear in p. As in the Euclidean case, collecting the coefficients 
of each pf, one observes that at a critical point ip is a coherent angle system. 
Also, if (f is a coherent angle system, then S h H yp does not depend on p. (Remem- 
ber that in the hyperbolic case, coherent angle systems satisfy ip e + tp- e < 8* 
instead of tp e + ip- e = 8* .) Again, this leads to a different variational principle. 



Proposition 10. If (p is a coherent angle system, then S^ yp = Sj^ p , wh 



ere 



+ Cl 2 (8* e +^ e + y_ e ) + Cl 2 (8* e -<p e - y- e ) - 2C1 2 (20|)). 

The sum is taken over non-oriented edges, or more precisely, over one repre- 
sentative for each pair of oppositely oriented edges. 

If \p is a critical point of S^ p under variations in the space of coherent angle 
systems, then the equations 

j sin ( g - ffff- ) .in = V; + V-' ) 
" = » ' 0S s ,„(»- + y . + .-.) s ,„(^ r ^)' m 

where f is the face to the left of e, coherently define p G M. F . This p is a critical 
point of S h yp. 

Conversely, suppose p is a critical point ofSh yp . Then equations (Q) define a 
coherent angle system, which is a critical point of S^ uc under variations in the 
space of coherent angle systems. 
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The proof is analogous to the proof of proposition ||. We would like to remark 
only on two points. 

First, equations ( |36| ) are equivalent to equations (||). This follows by a 
straightforward calculation. 



Second, to show that equations (36) are consistent if ip is a critical point of 
Sj^yp under variations in the space of coherent angle systems, proceed as follows. 
Suppose / is the face to the left of e, and e' = ere is the next edge in the 
boundary of /. It suffices to show that the value for pj obtained from the 
equation involving ip± e is the same as the value obtained from the equation 
involving ip± e >. The tangential space to the space of coherent angle systems 
is spanned by variations with tp e = —1, ip e t = 1 for two oriented edges e and 
e' = ere, and <p = otherwise. But dSj^ p /d(p e is twice the right hand side of 
equation (|3^) . This implies the consistency of equations (^ 

6.5 Leibon's functional 



Like Rivin, Leibon [ LeiOl considers circle patterns with arbitrary prescribed 
intersection angles. The circles correspond to the faces of a triangulation. 
Whereas Rivin treats the Euclidean case, Leibon treats the hyperbolic case. 
We believe that his functional can be identified with S^ yp . 

Appendix: The dilogarithm function and Clausen's 
integral 

In this appendix, we collect everything about the dilogarithm and Clausen's 
integral that is relevant for this paper. A more thorough treatment and an 



extensive bibliography are contained in Lewin's monograph [Lew81 



For \z\ < 1, the dilogarithm function is defined by the power series 
U 2 (z) 



z z 2 z 3 

12 + 22 + 32 



For \z\ < 1, 



and hence 



log(l - Z ) = - + — 



Li2(2) = _/'!2i(W), ( . (3r , 

Jo s 

In the light of this integral representation, one sees that the dilogarithm 
can be analytically continued to the complex plane cut from 1 to 00 along the 
positive real axis. 

Clausen's integral C^ix) can be defined by the imaginary part of the dilog- 
arithm on the unit circle: 

Cl 2 (z) =lmU 2 {e ix ) 

= i(Li 2 (e-)-Li 2 (e— )). 
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Clausen used the integral representation we will derive below to define this func- 
tion, hence the name. We consider Clausen's integral as a real valued function 
of a real variable. It is 27r-periodic and odd. The power series representation of 
the dilogarithm yields the Fourier series representation for Clausen's integral, 



Ch(x)= £ 



sin(na;) 



n=l 



Substitute C = ^ in the integral representation of the dilogarithm (^) to 
obtain, for < x < 2tt, 

Ch(x) = - jf log(2sin0 d£. 
Clausen's integral is almost the same as Milnor's Lobachevski function 
Jl{x) 



■J log|2sin£|e£=~Cl 2 (2z) 



Another related function appearing in the literature, also called Lobachevski 
function, is 



L(x) 



r i 

/ log(cos £) d£ = - Cl 2 (2x - tt) + x log 2 . 
Jo 2 



We will now derive a formula expressing the imaginary part of the diloga- 
rithm in terms of Clausen's integral not only on the unit circle, but anywhere 
in the complex plane. Suppose x is real and < 6 < 2ir. Substitute £ = + ^ 
in (|37|), to obtain 

ImL^e* + i6 ) = j % {^2{e x + l0 ) - Li 2 (e x ~ i6 )) 



1 



l-et- i0 \ (38) 



Now substitute 



and note that inversely 



Finally, one obtains 



1, fl-et- i6 
■log 



2z & ii_ e e 



e = io g 



2 sinr; 



2 sin(?7 + 9) 



Im Li 2 (e x + id )=yx+^ Cl 2 (2y) - ^ Cl 2 (2y + 20) + ± Cl 2 (20) , (39) 
where 

1 



2i l0g I 1 _ e x + i 
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From this, we derive the formula 



ImLi 2 (e a; + ^)+ImLi 2 (e- a; + ^) = 

px + Ch(p + d*) + Ch(-P + 0*) - Cl 2 (2e*) , (40) 

where 9* = tt — 9, and 

= _1_ (l + e x + i0 *)(l + e- x - i0 *) 
P ~ 2^ ° g (l + e ?-W*)(l + e- x + i0 *) ' 

Finally, note that p and x are related by 

tan(|)=tanh(|)tan( y ). 
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